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PREFACE 


In  late  1980,  a  study  to  produce  a  wave  climate  for  the  Mississippi 
Sound  region  was  initiated  at  the  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES).  This  study,  the  Mississippi  Sound  Wave-Hindcast 
Study,  was  authorized  by  the  U.  S.  Army  Engineer  District,  Mobile,  as 
part  of  a  larger  study  investigating  the  Mississippi  Sound  waters.  The 
Mobile  District  authorized  the  funds  throughout  the  study. 

The  report  describes  the  methodology  of  the  wind-field  development 
and  the  shallow-water  modeling  technique  and  presents  wave  hindcast  sta¬ 
tistics  for  the  period  of  1956-1975.  This  work  was  done  in  the  Hydrau¬ 
lics  Laboratory  under  the  direction  of  Mr.  H.  B.  Simmons,  Chief  of  the 
Hydraulics  Laboratory,  Dr.  R.  W.  Whalin,  former  Chief  of  the  Wave  Dy¬ 
namics  Division,  and  Mr.  C.  E.  Chatham,  Acting  Chief  of  the  Wave  Dynamics 
Division.  This  report  was  prepared  by  Dr.  R.  E.  Jensen.  Rebecca  Brooks, 
W.  D.  Corson,  D.  S.  Ragsdale,  and  H.  Messing  provided  additional  as¬ 
sistance  in  the  preparation  of  this  report. 

Commanders  and  Directors  of  WES  during  the  conduct  of  the  study 
and  the  preparation  and  publication  of  this  report  were  COL  Nelson  P. 
Conover,  CE,  and  COL  Tilford  C.  Creel,  CE.  Technical  Director  was 
Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units 

verted  to  metric  (SI) 

Multiply 

of  measurement  used 

units  as  follows: 

By 

in  this  report  can  be  con- 

To  Obtain 

feet 

0.3048 

metres 

knots  (international) 

0.514444 

metres  per  second 

miles  (U.  S.  nautical) 

1.852 

kilometres 

miles  (U.  S.  statute) 

1.609344 

kilometres 

miles  per  hour 

1.609344 

kilometres  per  hour 

miles  per  hour 

0.8689 

knots  per  hour 

square  feet  per  second 

0.09290304 

square  metres  per  second 

3 


MISSISSIPPI  SOUND  WAVE-HINDCAST  STUDY 


PART  I :  INTRODUCTION 

1.  The  primary  purpose  of  this  study  is  to  provide  a  20-year  wave 
climate  in  the  Mississippi  Sound  region  (Figure  1).  Secondary  efforts 
include  the  development  of  a  representative  wind  field  over  the  study 
area  and  the  development  of  a  wave  modeling  technique  that  provides  an 
accurate  representation  of  the  mechanisms  involved  in  shallow-water  wave 
growth  and  transformations. 

2.  The  sparseness  of  wave  information  in  Mississippi  Sound  be¬ 
comes  critical  to  the  design  of  coastal  or  offshore  structures  and  also 
to  the  operation  of  dredging  and  shipping  activities.  A  search  for 
available  wave  information  in  the  Sound  was  performed,  where,  over  a 
20-year  duration,  only  27  recorded  visual  ship  observations  reflected 
wave  conditions.  Wave  gaging  was  performed  with  limited  success  in  1980 
by  the  Raytheon  Corporation  (U.  S.  Army  Engineer  District,  Mobile,  1981) 
as  part  of  data  collection  effort  to  calibrate  and  verify  a  numerical 
model . 

3.  To  accurately  describe  the  wave  climate  within  Mississippi 
Sound,  two  types  of  techniques  can  be  employed--wave  gaging  and  wave 
hindcasting.  Although  a  network  of  wave  gages  might  eventually  provide 
a  good  data  source,  the  expense  involved  would  make  it  economically  pro¬ 
hibitive  to  provide  detailed  coverage  of  the  entire  Mississippi  Sound 
region.  Even  if  such  a  wave  gage  network  were  established,  there  would 
be  a  lag  time  of  many  years  before  sufficient  wave  information  was 
available  for  design  purposes. 

4.  A  viable  alternative  to  wave-gaging  Mississippi  Sound  would  be 
to  hindcast  the  wave  climate.  Many  different  techniques  are  presently 
available  to  generate  wave  characteristics.  These  techniques  can  be 
subdivided  into  three  categories:  empirical,  parametric,  and  numerical. 
All  three  types  of  wave-hindcasting  methods  have  certain  advantages  and 
disadvantages  associated  with  their  operation.  A  discussion  of  these 
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techniques  will  be  presented  later  in  this  report.  A  parametric  shallow- 
water  wave  model  was  developed  to  perform  the  20-year  hindcast  for  the 
Mississippi  Sound  region  to  minimize  computational  costs  without  loss 
in  describing  the  physics  of  the  problem. 

5.  The  use  of  hindcast  techniques  to  estimate  wave  heights  assumes 
that  the  coupling  between  the  atmospheric  boundary  layer  and  waves  gen¬ 
erated  by  the  motion  of  this  boundary  is  known.  Other  factors,  such  as 
dispersion  effects  (finite  water  depths),  bottom  effects  (friction),  and 
interactions  between  different  spectral  wave  components  must  also  be  con¬ 
sidered.  If  all  the  linear  and  nonlinear  wave  transformation  mechanisms 
are  adequately  described  and  if  the  atmospheric  motion  over  the  Sound  is 
known,  then  a  reliable  estimate  of  wave  conditions  can  be  made.  In  any 
estimate  of  wave  characteristics  (height,  period,  and  direction  of  wave 
propagation)  there  will  always  be  some  degree  of  error.  Because  of  the 
paucity  and  questionable  reliability  of  measured  wave  data  within  the 
Sound  (paragraphs  47-51),  an  extensive  error  analysis  cannot  be  made. 

All  assumptions  and  mechanisms  governing  the  hindcast  technique  have 
been  analyzed  in  the  "Wave  Data  Acquisition  and  Hindcast  for  Saginaw  Bay 
Study"  (Garcia  and  Jensen,  in  preparation).  Because  of  the  problems  as¬ 
sociated  with  the  measured  data  within  Mississippi  Sound,  the  error 
analysis  performed  in  the  Saginaw  Bay  study  is  presented  to  ensure  that 
the  modeling  technique  adequately  describes  the  processes  involved  in 
finite  water  depth  wave  generation. 

6.  PART  II  of  this  report  deals  with  the  wind  information  analy¬ 
sis,  PART  III  discusses  the  shallow-water  wave  modeling  technique, 

PART  IV  discusses  the  20-year  wave  hindcast  products,  and  PART  V  dis¬ 
cusses  results. 
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7.  The  recovery  and  subsequent  analysis  of  wind  data  become  an 
integral  part  in  the  Mississippi  Sound  Wave-Hindcast  Study.  The  lack  of 
continuous  long-term  wind  observations  over  the  region  forces  the  con¬ 
clusion  that  winds  over  Mississippi  Sound  must  be  estimated  from  alter¬ 
nate  sources  of  data.  Two  sources  of  data  with  a  sufficient  length  of 
record  (20  years)  are  synoptic  weather  maps  and  wind  observations  at  ad¬ 
jacent  land  stations.  The  conversion  of  a  synoptic  weather  map  from  its 
initial  state  to  a  finalized  wind  field  would  require  extensive  hand 
digitization  accompanied  by  an  inordinate  amount  of  computational  time. 
Therefore,  the  only  feasible  sources  of  data  for  the  estimation  of  winds 
over  Mississippi  Sound  are  from  land-situated  wind  observations  from 
data  sources  around  the  study  region. 

8.  Three  long-term  meteorological  sites  were  selected:  Keesler 
Air  Force  Base  (KAFB) ,  Biloxi,  Mississippi;  Mobile,  Alabama;  and  New 
Orleans,  Louisiana  (Figure  1).  The  primary  source  of  the  wind  informa¬ 
tion  was  taken  from  the  Biloxi  stations.  The  Biloxi  station  is  geo¬ 
graphically  the  closest  meteorological  station  to  the  study  area.  The 
weather  conditions  over  the  Sound  would  be  better  represented  employing 
this  data  set  over  the  wind  data  obtained  from  either  Mobile,  Alabama, 
or  New  Orleans,  Louisiana.  Since  gaps  in  the  wind  record  occurred,  it 
was  necessary  to  supplement  the  KAFB,  Biloxi,  data  with  either  the  Mobile 
or  New  Orleans  wind  information.  The  data  supplementing  procedure  im¬ 
plies  that  winds  in  the  area  are  nearly  uniform  both  in  intensity  and 
direction.  This  assumption  is  verified  in  Table  1  where  wind  conditions 
(speed  and  direction)  were  randomly  selected  (24  hr  of  data  per  set  per 
year).  If  the  KAFB,  Biloxi,  wind  data  are  assumed  to  be  the  norm,  the 
data  reveal  that  average  wind  speed  differences  of  12  percent  and  4  per¬ 
cent  are  found  at  Mobile  and  New_ Orleans,  respectively.  The  mean  wind 
direction  differences  are  approximately  3  to  4  deg  while  the  standard 
deviation  around  the  mean  is  36  to  44  deg.  The  large  variation  in  the 
standard  deviation  spans  approximately  two  wind  angle  classes  (winds  are 
given  in  16  direction  bands  at  22  5-deg  increments)  and  are  primarily 
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Table  1 

Differences  in  Wind  Speed  and  Direction! 


Location 
New  Orleans 
Mobile 


Wind  Speed 
Dif ferencett 
percent 

+4 

+  12 


Wind  Direction  Difference 


Mean 

Standard  Deviation 

deg 

deg 

-4 

±44 

-3 

±36 

t  Differences  between  cited  location  and  Biloxi, 
tt  All  wind  speeds  adjusted  to  33  ft. 


caused  by  low  wind  speed  conditions  wherein  the  wind  direction  becomes 
highly  variable  and  is  difficult  to  determine. 

9.  During  the  period  1956-1965,  wind  data  are  recorded  every 
hour.  The  wind  data  over  this  time  interval  are  averaged  over  a  3-hr 
interval  (using  the  wind  data  from  the  three  preceding  hourly  records). 
The  wind  data  from  1965-1975  are  given  in  3-hr  intervals;  thus  the  aver¬ 
aging  step  was  not  performed.  A  20-year  record  of  wind  speed  and  direc¬ 
tion  recorded  every  3  hr  is  generated  and  stored  on  magnetic  tape  for 
further  analysis. 

10.  Wind  speed  and  direction  changes  in  the  data  are  primarily 
caused  by  the  slow  migration  of  weather  patterns.  The  passage  of  frontal 
weather  patterns,  sea  breezes,  and  other  small-scale  meteorological  con¬ 
ditions  are  included  in  the  wind  data  although  their  overall  influence 

of  larger  scale  meteorological  conditions  would  be  diminished.  Also, 
tropical  storm  conditions  cannot  be  classified  within  this  constraint. 
During  the  20-year  hindcast  study  (1956-1975),  numerous  tropical  storms 
passed  through  or  passed  in  close  proximity  to  the  Mississippi  Sound 
region.  It  was  decided*  that  the  winds  would  be  treated  as  if  they  were 
generated  from  extratropical  storm  conditions.  These  data  are  shown  in 
Table  2  (Neumann  et  al.  1978). 

11.  Once  the  20-year  wind  data  are  assembled,  adjustments  must  be 


*  R.  Champion.  1982.  Personal  communication,  U.  S.  Army  Engineer  Dis¬ 
trict,  Mobile,  CE,  Mobile,  Alabama. 


Table  2 

Tropical  Storm  Condition  Data 


New 

Orleans 

Biloxi 

Mobile 

Max 

Max 

Max 

Wind 

Date 

Wind 

Date 

Wind 

Date 

Year 

Month 

Day 

knots 

day/hour 

knots 

day/hour 

knots 

day/hour 

1956 

06 

13-14 

19 

13/14 

23 

13/15 

22 

13/14 

1956 

09 

23-25 

42 

24/2 

31 

24/9 

35 

24/14 

1957 

09 

7-8 

28 

8/0 

12 

7/16 

17 

8/9 

1957 

09 

16-19 

22 

18/3 

29 

18/10 

30 

18/10 

1959 

05 

30-31 

24 

31/5 

19 

31/10 

20 

31/8-15 

1959 

10 

6-8 

14 

6/12 

14 

6/12 

18 

7/16 

1960 

09 

14-15 

15 

15/4-11 

36 

15/14 

30 

15/9 

1960 

09 

25-26 

14 

25/19 

14 

25/19 

20 

25/18 

1964 

09 

10-12 

14 

12/15 

13 

12/9 

20 

12/4 

1964 

10 

3-5 

35 

4/2 

33 

4/17 

35 

4/19 

1964 

11 

5-6 

17 

5/3-5 

5 

-- 

23 

5/12-15 

1965 

06 

14-16 

11 

15/15 

13 

14/13 

15 

14/15 

1965 

09 

8-12 

60 

10/0 

30 

9/23-10/1 

28 

10/0 

1965 

09 

28-29 

17 

29/0 

23 

29/14 

18 

29/18 

1969 

08 

16-18 

30 

17/21 

70 

17/23 

34 

18/0 

1969 

09 

5-6 

13 

5/15 

14 

5/15 

10 

6/9-12 

1969 

10-11 

30-1 

15 

30/21-1/3 

14 

1/4 

15 

1/3-6 

1970 

07 

21-22 

8 

-- 

10 

21/13 

12 

21/12 

1971 

09 

4-6 

14 

6/12 

* 

JL. 

13 

5/12 

1971 

09 

16-17 

25 

16/9-12 

* 

* 

24 

16/18 

1972 

06 

18-20 

17 

19/15 

* 

* 

16 

19/9-18 

1974 

09 

7-8 

27 

7/21 

15 

7/15-18 

19 

7/12 

1975 

09 

22-24 

23 

22/18 

30 

23/2 

21 

23/0 

*  Data  not  available. 
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made,  such  as  transforming  wind  speeds  from  the  given  elevation  to  the 
standardized  10-m  elevation,  converting  from  overland  to  overwater  wind 
speeds,  adjusting  for  air-sea  temperature  differences,  and  nonconstant 
drag  coefficients.  The  procedure  for  converting  given  wind  speeds  to 
"wave-model-ready"  winds  follows  that  described  in  Resio  and  Vincent 
(1976)  and  U.  S.  Army  Coastal  Engineering  Research  Center  (CERC), 

CETN-I-5  (1981). 

12.  No  long-term  sea  temperature  data  were  available  to  compute 
stable,  neutral,  or  unstable  air-sea  temperature  conditions.  It  was 
assumed  that  during  the  20-year  hindcast,  air-sea  conditions  were  un¬ 
stable  (as  described  in  U.  S.  Army  CERC  CETN-I-5),  thus  an  increase  of 
10  percent  was  added  to  the  resulting  wind  speed. 

13.  It  also 'must  be  noted  that  only  wind  conditions  coming  from 
the  west,  clockwise  to  the  east  (270  deg  <  0^  <  90  deg,  where  0^  is 
the  predominant  wind  direction),  were  converted  from  overland  to  ove*-- 
water  winds.  All  other  wind  conditions  come  from  overwater  areas 
relative  to  KAfB,  Biloxi,  as  shown  in  Figure  1,  and  are  not  adjusted. 

14.  In  summary,  the  20-year  wind  conditions  are  derived  from  a 
single  source  (Biloxi,  Mississippi,  KAFB),  and  supplemented  with  two  ad¬ 
ditional  land  stations  (Mobile  and  New  Orleans).  The  winds  over  the 
study  area  were  shown  to  be  uniform  in  speed  and  direction  (Table  1). 

Wind  conditions  are  then  converted  from  overland  to  overwater  winds  and 
adjusted  for  a  nonconstant  drag  coefficient  and  unstable  air-sea  tempera¬ 
ture  differences.  From  this  analysis  a  continuous  20-year  record  is  ob¬ 
tained  of  surface  winds  derived  from  all  available  measured  wind  informa¬ 
tion  in  the  Mississippi  Sound  region. 
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PART  III:  SHALLOW-WATER  WAVE  MODELING  TECHNIQUE 


Introduction 

15.  The  predictions  of  shallow-water  wave  characteristics  have 
become  a  focal  point  of  research  activities  across  the  world.  Because 
construction,  shipping,  and  dredging  operation  costs  have  drastically 
increased  over  the  years,  coastal  engineers  have  been  faced  with  more 
accurately  defining  the  shallow-water  wave  climate.  A  better  under¬ 
standing  of  shallow-water  wave  growth  and  transformation  mechanisms  is 
slowly  evolving  through  controlled  wave-measuring  programs  such  as  the 
Atlantic  Remote  Sensing  Land-Ocean  Experiment  (Vincent  and  Lichy  1981). 
However,  not  all  of  the  questions  have  been  answered,  and  it  will  take 
some  time  before  all  shallow-water  wave  transformation  mechanisms  are 
quantified.  In  light  of  this,  the  shallow-water  wave  modeling  technique 
(SWWM)  employed  in  this  study  adopts  "state-of-the-art"  mechanisms  cur¬ 
rently  available.  The  main  intent  in  the  development  of  the  SWWM  is  to 
describe  the  physical  processes  as  accurately  as  possible  while  simpli¬ 
fying  the  computational  procedures  to  a  degree  where  a  long-term  hindcast 
study  is  economically  feasible. 

16.  This  part  of  the  report  is  subdivided  into  three  sections: 

(a)  the  discussion  of  the  theoretical  aspects  of  the  SWWM,  (b)  the  model 
setup  applied  specifically  to  Mississippi  Sound,  and  (c)  calibration  and 
verification  to  existing  measured  and  observed  wave  conditions. 

17.  The  existing  shallow-water*  wave  prediction  techniques  can  be 
classified  into  three  categories--empirical ,  numerical,  and  parametric. 

A  brief  review  is  given  below,  but  a  more  detailed  description  of  these 
techniques  can  be  found  in  Hsiao  (1978)  and  Vincent  (1982). 

18.  The  most  widely  used  technique  to  generate  shallow-water  wave 
conditions  is  the  empirical  method  first  described  by  Bretschneider 
(1954).  Since  that  time,  this  approach  has  been  revised  (CERC  1977, 

*  Shallow  water  is  defined  as  water  depth  conditions  where  dispersion 
effects  become  important  (or  changes  in  the  phase  and  group  wave 
celerities  become  dependent  on  changes  in  the  water  depth). 


CERC  CETN-I-6  1981).  The  significant  wave  height  (H  )*  and  significant 
wave  period  (T  )  are  estimated  by  empirical  relationships  between  wave 
height  and  wind  speed,  duration,  and  fetch.  Also  included  in  this 
modeling  technique  is  the  influence  of  bottom  friction  mechanisms  on 
wave  dissipation  formulated  by  Bretschneider  and  Reid  (1954). 

19.  The  Bretschneider  empirical  shallow-water  wave  method  is 
based  only  on  experimental  data  and  no  consideration  of  the  physical 
processes  involved  in  wave  development,  such  as  generation  and  interac¬ 
tions,  is  included  in  the  technique. 

20.  Collins  (1972)  developed  a  finite  depth  wave  model  by 
numerically  solving  the  energy  balance  equation.  The  model  included  the 
effect  of  wind  generation  similar  to  that  employed  by  Barnett  (1968). 
Also,  the  transformation  mechanisms  modeled  included  wave  refraction, 
shoaling,  bottom  friction,  and  wave  breaking.  Another  shallow-water 
wave  model  based  on  the  ray  technique  was  developed  by  Cavalri  and 
Rizzoli  (1981).  Their  model  included  the  same  processes  as  Collins' 
version.  One  process  that  was  not  considered  in  either  model  was  the 
incorporation  of  the  nonlinear  energy  transfers,  or  wave-wave  interac¬ 
tions  (Hasse lmann  1962).  The  influence  of  this  mechanism  has  been  shown 
to  bf  amplified  in  a  finite  water  depth  wave  environment  (Heterich  and 
Hasselmann  1980). 

21.  A  parametric  wave  prediction  model  was  first  developed  by 
Hasselmann  et  al.  (1976)  for  deepwater  wave  conditions.  More  recently 
Hsiao  (1978)  and  Shemdin  et  al.  (1980)  have  adopted  the  original  work 
of  Hasselmann  et  al.  (1976)  to  finite  water  depths.  The  effect  due  to 
wind  generation,  refraction,  shoaling,  nonlinear  bottom  interactions 
(friction,  percolation,  viscoelasticity,  and  scattering),  and  nonlinear 
wave-wave  interactions  (Heterich  and  Hasselmann  1980)  were  formulated  in 
the  model.  The  inhomogeneous  wave  equation  was  solved  numerically  in 
one  dimension  using  a  Runge-Kutta  (numerical  approximation  method)  type 
of  approach  of  successive  approximations. 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  de 
fined  in  the  Notation  (Appendix  B). 
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22.  No  direct  measurements  can  be  made  to  verify  the  relative  im¬ 
portance  of  each  transformation  mechanism;  thus  calibration  and  verifica¬ 
tion  of  this  complex  parametric  model  can  be  very  difficult.  One  has  to 
consider  parameters  at  each  selected  site  (such  as  bottom  topography  and 
sediment  size  and  type)  and  continually  make  adjustments  to  the  model 

for  changes  in  time  at  that  site.  Therefore  this  technique  is  restricted 
to  the  computation  of  site-specific  modeling  and  is  difficult  to  employ 
for  large  geographical  areas  unless  extensive  measured  data  exist. 

23.  Since  the  above  approaches  to  an  accurate  physical  estimation 
of  wave  conditions  in  a  large  geographical  area  are  inappropriate,  an 
alternate  approach  to  the  problem  of  long-term  wave  hindcasting  was  de¬ 
veloped  specifically  to  meet  the  needs  of  the  Mississippi  Sound  Wave- 
Hindcast  Study.  It  should  be  recognized  at  the  outset  of  this  discus¬ 
sion  that  in  order  to  maximize  computational  efficiency,  a  new  approach 
to  shallow-water  wave  modeling  was  developed.  The  key  to  this  new  ap¬ 
proach  is  that  the  resulting  wave  conditions  are  produced  by  transforma¬ 
tion  mechanisms  rather  than  transforming  wave  conditions  during  wave 
propagation.  Computationally,  it  is  more  efficient  to  transform  wave 
conditions  from  one  location  to  another  than  it  is  to  simultaneously 
propagate  and  transform  sets  of  wave  conditions  from  a  given  location  to 
the  next.  Thus  a  shallow-water  wave  modeling  technique  was  developed 
that  includes  wave  growth,  nonlinear  transfers  of  energy,  shoaling,  re¬ 
fraction,  energy  dissipation  resulting  from  wave  decay  (if  during  growth 
the  spectrum  is  fully  saturated),  bottom  friction,  and  wave  breaking. 

Most  transformation  mechanisms  are  evaluated  parametrically,  while  the 
remaining  mechanisms  are  determined  from  empirically  derived 
relationships . 

24.  Prior  to  the  discussion  of  the  transformation  mechanisms,  it 
is  necessary  to  discuss  the  assumptions  governing  the  SWWM  results.  Wind 
conditions  are  assumed  to  be  uniform  over  Mississippi  Sound.  These  con¬ 
ditions  are  typical  of  those  generated  from  slowly  moving  extratropical 
storm  systems.  The  winds  also  are  assumed  to  remain  constant  for  each 
3-hr  observation.  Wave  propagation  is  assumed  to  be  restricted  to  the 


13 


direction  described  by  the  winds.  Over  each  given  fetch  length  (de¬ 
pendent  on  a  wind  angle  class),  the  bottom  topography  is  represented  by 
straight  and  parallel  bottom  contours.  Mississippi  Sound  is  assumed  to 
be  a  closed  system,  whereby  all  wave  conditions  found  in  the  Sound  are 
generated  within  its  boundaries.  No  energy  is  permitted  to  propagate 
through  the  barrier  island  system  to  the  south  or  from  the  east,  due  to 
Mobile  Bay  (Figure  1).  Without  a  comprehensive  wave  gaging  study  in  and 
around  the  Mississippi  Sound  region,  an  accurate  assessment  of  this  as¬ 
sumption  cannot  be  made.  The  intent  of  the  given  assumptions  is  not  to 
overly  simplify  the  generation  of  a  long-term  hindcast  for  Mississippi 
Sound,  but  make  it  economically  feasible  without  restricting  the  overall 
physics  of  the  problem. 


Theory 

25.  This  section  describes  the  theory  employed  in  the  SWWM. 

There  are  five  basic  parts  to  the  SWWM:  (a)  wave  growth,  (b)  applica¬ 
tion  of  conservative  and  nonconservative  mechanisms  on  the  growing 
energy,  (c)  construction  of  a  one-dimensional  energy  density  spectrum 
based  on  available  energy,  (d)  application  of  an  additional  atmospheric 
source  term  on  the  spectrum,  and  (e)  limitation  of  the  energy  due  to 
wave  breaking. 

26.  Hasselmann  et  al.  (1976)  introduced  a  parametric  model  of 
wind-wave  generation,  relating  the  rate  of  energy  growth  to  nondimen- 
sional  characteristics  of  the  wind  field.  This  energy  growth  (in  space 
or  time)  is  governed  by  a  self-similar  process  and  verified  through  ex¬ 
tensive  prototype  data  (Hasselmann  et  al.  1973  and  1976).  In  these 
studies,  the  dominant  energy  input  on  the  forward  face  of  the  spectrum 

is  related  to  convergence  of  energy  flux  due  to  nonlinear,  resonant  wave- 
wave  interactions  (Figure  2)  of  the  form  described  by  Hasselmann  (1962). 
Work  conducted  by  Mitsuyasu  (1968,  1969)  and  Kitaigordskii  (1962)  also 
displayed  similar  results.  Although  these  formulations  are  governed  for 
deepwater  wave  conditions,  they  are  used  in  the  SWWM  simply  because  the 
only  formulation  of  the  nonlinear  transfers  in  shallow  water  are  based 
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Figure  2.  Schematic  representation  of  the  nonlinear 
wave-wave  interaction 


specifically  on  JONSWAP-type  wave  spectra.  The  parameterizations  of  the 
shallow-water  nonlinear  transfers  become  too  complex  and  require  addi¬ 
tional  computational  time  to  resolve.  Also,  the  JONSWAP-type  spectra 
contain  three  free  parameters  which  are  based  primarily  on  open-ocean, 
deepwater  spectra  wave  conditions  that  are  not  appropriate  to  Mississippi 
Sound. 

27.  The  rate  of  wave  growth  under  ideal  conditions  of  fetch  limi¬ 
tation  or  duration  limitation,  governed  by  a  stationary  wind  field,  can 
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be  computed  from  Hasselmann  et  al.  (1976).  For  growth  along  a  fetch, 
the  solution  yields 


-7  2  F 

E  =  1.6  x  10  U  - 
o  g 


and  for  growth  through  time,  it  becomes 


(1) 


E 

o 


4.3  x  10 


-10^8/7-4/7^0/7 


(2) 


where 

E  =  total  energy  resulting  from  a  wind  speed  U  (assumed  to  be 
overwater  wind  conditions  adjusted  to  33  ft*  in  elevation), 
blowing  over  a  given  fetch  length  F 

g  =  acceleration  due  to  gravity 

t  =  total  elapsed  time  since  the  wind  began  to  blow 
28.  Two  additional  pieces  of  information  are  required  to  quantify 
the  distribution  of  Eq  given  in  the  form  of  an  energy  density  spectrum. 
The  nondimensional  peak  frequency,  f  ,  and  the  Phillips’  equilibrium 
constant,  a  (Phillips  1957),  are  shown  in  Figures  3  and  4.  Mathe¬ 
matically  these  parameters  are  represented  by 


and 


— 0  22 
a  =  0.076X 


(3) 


f 

m 


~  —  0.33 
3.5X 


where  X  is  the  nondimensional  fetch  length 


(4) 


(5) 


29.  The  selection  of  either  fetch  (Equation  1)  or  duration 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ments  to  metric  (SI)  units  is  presented  on  page  3. 
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Figure  3.  Nondimensional  peak  frequency  as  a 

function  of  nondimensional  fetch  length  X  (after 
Hasselmann  et  al.  (1973)) 


XT'  i  io  »'  to*  to*  to* 

PtMftip*  constant  v.  fetch  scaled  according  to  KAaigordskn.  SmoK- fetch  data  art  ototamed  from  wind-ware  tonka. 
(CapMary-wave  data  waa  exchidod  whore  poeaiblo.)  Maaauramante  by  Sutherland  (1967)  and  Tobo  (1971)  wore 
taken  from  W  J.  Pi  arson  and  RA.  Stacy  (19731 

Figure  4.  Phillips'  equilibrium  constant  or  as  a  function  of 
nondimensional  fetch  length  X  ,  after  Hasselman  et  al.  (1973) 


limited  conditions  is  determined  from  the  following: 


t  . 
mm 


(6) 


where 


t  .  =  minimum  duration  condition 

min 

Tg  =  significant  wave  period  (CERC  CETN-I-6  1981) 
given  by: 


T  = 
s 


(7) 


where  h  is  the  mean  water  depth  along  F  . 

30.  If  t  .  is  less  than  3  hr  (duration  of  each  input  wind  con- 

min  r 

dition)  then  Equation  2  will  be  used  to  compute  the  total  energy;  other¬ 
wise,  Equation  1  will  be  employed. 

31.  The  parameterization  of  the  wave  growth  is  somewhat  restricted 
such  that  when  the  nondimensional  peak  frequency  attains  a  value  of  0.13 
or  less,  a  fully  developed  sea  state  is  achieved  and  wave  growth  is 
automatically  halted.  Over  long  fetch  lengths  and  low  wind  speeds  this 
condition  can  occur  to  some  degree  of  regularity.  Thus  Equations  1-5 
are  then  redefined  by 


10 

Q  =  *  I 

i=l 


(8) 


where 

Q  =  the  dependent  parameters 
K  =  the  nonvarying  parameters  (and  constants) 

=  the  independent  parameters  (F  and  X)  found  in  Equations  1-5 
The  parameter  i  is  the  increment  counter.  After  each  discrete  fetch 
length  F .  ,  the  nondimensional  peak  frequency  is  evaluated  to  determine 
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if  f  <  0.13  .  If  this  occurs,  wave  growth  is  curtailed,  and  wave  de¬ 
ni  — 

cay  is  initiated  for  the  remainder  of  the  fetch  length.  Wave  decay  is 

parameterized  following  the  work  conducted  by  Bretschneider  (1952)  and 

Mitsuyasu  and  Kimura  (1965)  for  f  (where  f  =  f  g/U)  while  the  total 

m  m  m 

energy  decay  rate  follows  that  described  by  Jensen  (in  preparation). 

32.  Wave  conditions  generated  in  Mississippi  Sound  must  also  con¬ 
sider  dispersion  effects  resulting  from  finite  water  depth  conditions. 

When  the  water  depths  vary  from  F.  to  F.  ,  ,  one  must  consider  the 

i  l+l  ’ 

conservative  transformation  mechanisms  of  shoaling  and  refraction.  Wave 
shoaling  is  determined  from  the  evaluation  of  group  celerities  governed 
by  linear  theory.  Wave  refraction  is  neglected  under  the  assumptions 
that:  (a)  waves  travel  in  the  direction  of  the  winds,  and  more  impor¬ 

tantly  (b)  the  bottom  topography  is  assumed  to  be  straight  and  parallel 
(perpendicular  to  fetch  direction)  for  every  fetch  length  (discussed  in 
the  methodology  section).  The  latter  assumption  constrains  only  limited 
portions  of  the  Mississippi  Sound  region,  where  wave  propagation  can,  in 
its  extreme  case,  travel  parallel  to  the  bottom  contours.  Considering 
the  water  depths  in  Mississippi  Sound  and  peak  wave  periods  (T  =  1 / f m ) 
in  the  range  of  2  to  8  sec,  wave-refraction  effects  (and  subsequent 
"errors")  would  be  about  2  to  25  percent  as  shown  by  the  crosshatched 
area  in  Figure  5.  This  is  assuming,  at  most,  that  the  angle  between  the 
wave  crest  and  bottom  contour  is  30  deg.  The  initial  direction  of  wave 
propagation  is  limited  to  16  angle  classes  (because  of  the  wind  data  em¬ 
ployed  in  this  study);  thus  the  accuracy  in  the  resultant  refracted  wave 
condition,  by  similarity,  cannot  be  any  more  accurate  than  initially 
specified.  Any  differences  in  the  resulting  wave  climate  directional 
distribution  would  be  minimal. 

33.  Finite  water  depth  conditions  also  lead  to  bottom  dissipation 
effects  on  the  growing  seas.  Energy  losses  associated  with  bottom  fric¬ 
tion  are  empirically  modeled  using  the  following  sets  of  equations  de¬ 
veloped  by  Bretschneider  and  Reid  (1954): 


E 


ffE.^.AF.f4 
If  i  m 

K 


-1 


(9) 
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Figure  5.  Changes  in  wave  direction  due  to  refraction  on  slopes  with 
straight  and  parallel  depth  contours.  The  crosshatched  area  defines 
the  range  of  wave  conditions  typically  found  in  Mississippi  Sound 


where 

E  =  the  final  total  energy  at 
Ej  =  the  original  total  energy  at  F._j 

ff  =  the  dimensionless  friction  factor  (set  for  Mississippi  Sound 
at  0.001) 

AF.  =  the  distance  of  wave  travel  within  the  discrete  fetch 
length,  where 


and 


(10) 


♦f 


(ID 
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where 


kj  =  the  wave  number  (k^  =  2rt/L^) 

L.  =  the  wavelength  evaluated  for  f 
1  m 

h^  =  the  water  depth  at 

34.  The  second  theoretical  aspect  of  SWWM  deals  primarily  with  the 

distribution  of  the  total  energy  (E  )  in  the  form  of  a  one-dimensional 

o 

discrete  frequency  spectrum  E(f^)  •  Through  the  use  of  similarity 
principles,  Kitaigordskii ,  Kratsitskii,  and  Zaslavaskii  (1975)  extended 
Phillips'  deepwater  hypothesis  (Phillips  1958)  of  the  equilibrium  range 
in  the  spectrum  of  wind-generated  surface  waves  to  finite  depth  condi¬ 
tions.  The  spectral  form  is  defined  by 


E(f..)  =  ag2(2rc)”/*f”5<t>(mh)  f.  >  fm 


(12) 


where 


E(f.)  =  the  energy  density  at  each  discrete  frequency  band,  f. 
J  J 

g  =  acceleration  due  to  gravity 

=  a  nondimens ional  function  dependent  on  given  by 

J/2 


% 


=  2nf 


j(!) 


(13) 


The  function  'Ku^)  varies  from  1.0  in  deep  water  to  0.0  when  h  =  0.0  , 
as  shown  by  Figure  6.  When  is  less  than  1.0,  ♦  (tu^)  can  be  ap¬ 

proximated  by: 


u>. 


(14) 


and  therefore, 


E(f.)  =  \  agh(27i)"2fT3  f.  >  fm  (15) 

-5  -3 

or,  the  spectral  shape  changes  from  an  f  to  an  f  in  the  tail  of 
the  energy  density  spectrum  and,  more  importantly,  becomes  a  function  of 
the  water  depth. 
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Figure  6.  The  universal  dimensionless  function  <t>  (solid  curve) 
2 

and  the  function  ui^/2  (dashed  curve)  (after  Kitaigordski i , 
Krasitskii,  and  Zaslavaskii  (1975)) 


by 


35.  The  forward  face  of  the  spectrum  is  assumed  to  be  represented 


E(f.)  =  cxg2(2w)‘V5  exp 
j  m  r 

where  <t>'(u),  )  is  evaluated  from  the  w,  defined  at  f  .  Field  and 
h  h  m 

laboratory  data  by  Coda  (1974),  Thornton  (1977),  Ou  (1980),  Iwata  (1980), 
and  Vincent  (1981)  support  the  form  given  by  Equation  15.  The  verifica¬ 
tion  of  Equation  16  can  be  found  in  Vincent"  and  Garcia  and  Jensen  (in 
prepa  ration) . 

36.  The  parametric  representation  of  wave  growth  assumes  a 
dynamic  balance  between  atmospheric  sources  and  transfers  of  energy  re¬ 
sulting  from  wave-wave  interactions  (Figure  2).  This  parameterization 
was  based  on  deepwater  wave  conditions,  Hasselmann  et  al .  (1976).  Dur¬ 
ing  the  "Wave  Data  Acquisition  and  Hindcast  Study  for  Saginaw  Bay"  it 
was  determined  that  over  moderately  short  fetch  lengths  (10  to  20  nm) , 


1  - 


<J>' (wh) 


f.  <  f 
J  m 


(16) 


*  C.  L.  Vincent  1982a.  Personal  communication,  U .  S.  Army  CERC , 
Fort  Belvoir,  Va . 
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this  deepwater  growth  rate  expression  (Equations  1  and  2)  consistently 
underpredicted  the  total  energy  found  in  the  measured  data  (Garcia  and 
Jensen,  in  preparation).  The  only  theoretically  consistent  location  to 
add  the  energy  would  be  on  the  forward  face  of  the  spectrum  (Figure  7). 


The  function  E(f,h) 


THEORY 


is  the  saturated  spectrum  based  on  Equa¬ 


tions  12  and  16,  and  E(f,h) 


WEIGHTED 


is  the  spectrum  based  on  E 


after  wave  growth.  This  process  also  shifts  f  to  a  lower  frequency 
which  has  been  noticed  in  field  data  (Vincent*).  As  the  fetch  length 
increases,  the  relative  amount  of  added  energy  decreases;  and  eventually, 
no  additional  energy  is  incorporated  into  the  resulting  spectrum. 

37.  It  has  been  shown  that  the  water  depth  greatly  influences  the 
spectral  shape  and  in  so  doing  will  influence  the  maximum  wave  condition. 
The  parametric  formulation  follows  the  work  conducted  by  Vincent  (1981). 
The  depth-limiting  maximum  wave  condition  is  given  by, 


E(f )  df 


where 


H  =  maximum  wave  condition 
m 


f  =  lower  frequency  bounding  the  total  energy  (equal  to 

C  0.9f  ) 

m 

E(f)  =  from  Equation  12 

Integrating  Equation  17,  one  obtains  the  absolute  limit  on  the  wave  con¬ 
dition  at  a  particular  water  depth,  where 

,  .,1/2 

H  =  (SBl't—  (18) 

m  ni 

c 

38.  In  summary,  the  physical  process  governing  wave  generation  and 
transformations  has  been  theoretically  determined  using  available,  state- 
of-the-art  techniques.  It  must  be  emphasized  that  not  all  shallow-water 


*  C.  L.  Vincent  1982b.  Personal  communication,  U.  S.  Army  CERC, 
Fort  Belvoir,  Va. 
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Figure  7.  Construction  of  the  final  energy  density  spectrum 
(solid  symbols)  caused  by  shallow-water  wave  generation 


transformation  processes  have  been  (or  can  be)  measured  to  determine 
their  relative  effect  on  the  total  energy,  spectral  shape,  and  peak  fre¬ 
quency.  Therefore,  the  development  of  the  SWWM  employed  in  the  Missis¬ 
sippi  Sound  Wave-Hindcast  Study  attempts  to  model  the  physics  of  the 
problem  in  a  general  sense,  while  effectively  maximizing  computational 
efficiency. 


SWWM  Setup  for  Mississippi  Sound 

39.  The  SWWM  is  designed  to  compute  wave’  conditions  at  selected 
locations  in  Mississippi  Sound.  All  wave  conditions  generated  assume 
constant  water  depths  over  time,  therefore  neglecting  changes  in  the 
water  elevation  caused  by  tides,  freshwater  discharges,  and  surges.  In 
order  to  improve  computational  efficiency,  a  polar  coordinate  system  is 
selected  wherein  the  origin  is  placed  at  each  of  the  23  selected  station 
locations  (Table  3  and  Figure  8).  Fetch  length  rays  are  projected 
outward  from  the  origin  at  20-deg  intervals.  A  total  of  18  rays  exist 
for  each  station.  The  selection  of  the  18  rays  at  20-deg  intervals  as¬ 
sures  that  the  variability  of  the  shoreline  boundaries  (other  than  0,  90, 
180,  and  270  deg)  are  accurately  described  rather  than  using  only  the  16 
wind  angle  class  directions.  Fetch  lengths  and  water  depths  are  dis¬ 
cretized  into  10  subsections  along  the  total  length  of  every  ray.  The 
water  depth  selected  for  each  subsection  is  interpolated  from  available 
NOAA  bathmetric  charts.  The  parameters  h^  (discrete  water  depth),  F^ 

/  10  \ 

(discrete  fetch  length),  and  F  (total  fetch  length,  where  F  =  £  F.J 

1  \  i=l  7 

then  become  direct  functions  of  a  given  wind  direction.  This  is  shown 

schematically  for  Station  1  in  Figures  9  and  10.  All  fetch  lengths  and 
water  depths  are  tabulated  in  Appendix  C.  Special  fetch  lengths  are 
provided  for  certain  station  locations  when  the  given  wind  direction  is 
either  90-deg  or  270-deg  azimuth  (see  Appendix  C).  In  many  instances 
the  wind  direction  will  not  correspond  identically  to  a  given  fetch 
length  ray.  When  this  occurs,  a  new  ^(0^)  and  h^(F^,0w)  are  com¬ 
puted  via  linear  interpolation  between  two  discrete  fetch  length  rays. 
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Table  3 

Station  Locations  and  Water  Depths 


Station 

No. 

Longitude,  W 
deg 

Latitude,  N 
deg 

Water  Depth 
ft 

i 

88.17 

30.28 

6.5 

2 

88.25 

30.28 

9.5 

3 

88.30 

30.20 

12.0 

A 

88.30 

30.25 

16.5 

5 

88.  A2 

30.25 

17.0 

6 

88.50 

30.25 

19.0 

7 

88.58 

30.25 

10.0 

8 

88.58 

30.32 

8.0 

9 

88.67 

30.29 

13.0 

10 

88.75 

30.29 

11.5 

11 

88.83 

30.33 

9.0 

12 

88.83 

30.25 

15.0 

13 

88.92 

30.33 

11.0 

1A 

88.92 

30.25 

15.0 

15 

89.00 

30.33 

10.5 

16 

89.00 

30.25 

16.0 

17 

89.08 

30.33 

9.0 

18 

89.08 

30.26 

15.0 

19 

89.17 

30.29 

13.0 

20 

89.25 

30.25 

9.0 

21 

89.25 

30.17 

12.5 

22 

89.33 

30.21 

11.0 

23 

89.33 

30.15 

11.0 

26 


Schematic  representation  of  the  polar  coordinate  system  employed  In  the 
Mississippi  Sound  Wave-Hlndcas t  Study 


FETCH  INCFIEMENT 


Figure  10.  Schematic  representation  of  the  discrete  water  depths  as  a 

function  of  wind  direction 

AO.  For  the  shallow-water  wave  spectrum  E(f.)  defined  in  the 
previous  section,  35  frequency  bands  are  selected  that  will  adequately 
describe  the  energy  distribution.  The  frequency  bands  range  from  0.083 
to  2.0  Hz  and  are  tabulated  in  Appendix  C.  Since  (Equations  12 

and  16)  is  dependent  on  each  discrete  frequency  and  the  water  depth  at  a 
given  station  location,  all  values  were  then  precomputed  and  stored  on 
file  to  be  used  later  in  the  computational  scheme. 

41.  The  input  conditions  to  the  SWWM  are  the  wind  speed  (given  at 
33  ft  measured  overland)  and  the  wind  direction  (measured  "from-which- 
they-came,"  in  degrees  azimuth).  The  wind  speed  is  then  adjusted  to 
overwater  winds  as  previously  discussed  in  PART  II.  The  proper  fetch 
length,  F^  ,  and  water  depths,  h^  ,  are  then  selected  for  each  station 
for  the  given  wind  direction.  With  this  information,  the  SWWM  is  now 
ready  to  compute  the  wave  conditions  throughout  Mississippi  Sound.  This 
procedure  is  followed  for  every  3-hr  interval  over  the  20-year  period  of 
record. 

42.  With  all  dependent  parameters  known  (^0^,  *^(6^),  h^(F^,  9w)^ 


29 


wave  growth  can  be  initiated.  The  first  step  in  the  growth  sequence 
selects  either  duration  or  fetch-limited  conditions.  From  that  point, 
total  energy  growth  is  dependent  on  a  discrete  increment  of  time  (t^)  Jr 
fetch  length  (F^).  Wave  transformation  processes  then  act  only  on  the 
portion  of  energy  created  in  a  given  fetch  (or  time)  increment.  Changes 
in  the  group  speed  (h^,  f^)  are  determined  from  the  changes  in 

water  depth  from  F^  j1  to  F^  ,  and  assuming  f^  remains  constant. 
Energy  losses  associated  with  bottom  friction  effects  are  also  treated 
in  a  similar  manner.  Figure  11  schematized  the  processes  involved  in 
fetch-limited  energy  growth.  Duration-limited  energy  growth  follows  the 


same  procedure  where  temporal  growth  rates  are  employed.  When  satura¬ 
tion  conditions  exist  (f  <  0.13),  energy  growth  will  terminate  and  the 
energy  will  be  decayed.  During  the  latter  sequence,  the  transformation 


processes  will  continue  to  act  on  the  residual  energy 


(ERES^ 


following 


the  previously  described  procedure. 


43.  Once  E  is  determined,  it  now  becomes  a  matter  of  distrib- 
o 

uting  that  energy  over  E(fj)  •  This  is  performed  at  each  station,  em¬ 
ploying  a  slightly  modified  version  of  Equations  12-16.  Integrating 

E(f.,  h  )  ,  now  dependent  on  the  water  depth  at  a  given  station  (h  ), 

J  s  s 

one  finds  the  total  saturated  energy  caused  by  a  given  wind  condition 


blowing  over  a  fetch  length  for  some  period  of  time: 


E 

°SAT 


h  )  of 
s 


(19) 


Therefore,  the  resulting  spectrum  E*(f.,  h  )  is  then  scaled  according 

J  s 

to  Eq  (the  total  available  energy).  Hence, 

E* ( f  . ,  h  )  =  yE ( f  .  ,  h  )  (20) 

J  S  J  ® 

where  y  is  equal  to  E  /F. 

0  °SAT 

44.  The  resulting  spectrum  E*(fj,  hR)  is  then  readjusted  to 
compensate  for  employing  a  deepwater  parametric  form  of  energy  growth. 
The  procedure  follows  the  technique  described  earlier  in  this  section. 
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If  conditions  are  such  that  the  total  energy  is  decaying  during  wave 
growth,  this  step  is  neglected.  Thus,  at  this  point  in  the  SWWM,  a  one- 
dimensional  discrete  energy  density  spectrum  is  provided  for  each  of  the 
23  stations  for  one  3-hr  period. 

45.  The  wave  parameters  are  computed  by 


and 

T  =  jr-  (22) 

m 

where  f^  is  equal  to  f^  ,  where  the  maximum  energy  density  occurs. 

The  direction  of  wave  propagation  <)>  is  assumed  to  be  equal  to  the 
wind  direction  for  all  stations. 

46.  A  final  check  is  made  to  determine  if  H  (Equation  21)  is 

s 

greater  than  H  (Equation  18).  If  H  >  H  ,  then  E*(f.,  h  )  must 
&  m  ^  s  m  ’  j’s 

be  scaled  according  to  H  similar  to  Equation  20,  where  y  now  equals 
(»>0  • 

Calibration  and  Verification 

47.  In  all  wave  hindcasting  studies,  comparisons  to  actual  gage 
measurements  are  necessary  to  determine  the  validity  of  the  long-term 
wave  statistics.  It  is  unfortunate  that  ship  observations  and  wave  gage 
data  in  Mississippi  Sound  are  very  limited.  In  a  search  of  20  years  of 
ship  observation  records,  only  27  records  were  obtained.  Although  a 
wave  gaging  program  was  performed  by  Raytheon  Ocean  Systems  Company  as 
part  of  the  Mississippi  Sound  Data  Collection  Program,  only  one  of  the 
two  wave  data  sets  prepared  was  acceptable  for  comparisons  to  the  SWWM 
results . 

48.  The  SWWM  results  are  compared  with  the  limited  ship  records 
during  the  period  of  3-10  March  1965  (Figure  12).  It  should  be  noted 
that  measured  data  are  incremented  in  0.5-m  intervals  and  the  location 


WAVE  DIRECTION.  DEG  AZM 
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of  the  ship  is  available  only  to  the  nearest  0.1  deg.  Both  of  these 
factors  make  one-to-one  comparisons  difficult.  Overall,  the  SWWM 
results  display  similar  trends  observed  in  the  ships  records  for  both 
4*  and  H 

s 

49.  As  mentioned,  the  wave  gaging  program  involved  in  the  Missis¬ 
sippi  Sound  reconnaissance  study  (USAED,  Mobile,  1981)  had  limited  suc¬ 
cess.  Review  of  the  wave  gage  operations  revealed  that  the  gage  could 
not  resolve  wave  periods  less  than  5.0  sec.  Since  Roberts*  concluded 
that  61  to  71  percent  of  the  waves  south  of  Mississippi  Sound  have  wave 
periods  of  less  than  5.0  sec  and  since  the  gage  could  not  resolve  these 
conditions,  comparisons  of  hindcast  and  gage-measured  wave  periods  are 
not  appropriate  for  verification. 

50.  The  only  station  at  which  suitable  prototype  measurements  of 

wave  height  are  available  is  actually  outside  of  Mississippi  Sound  just 

south  of  Ship  Island  (Figure  8,  WG  4).  Comparisons  of  hindcast  and  gage- 

measured  H  are  presented  in  Figure  13.  The  hindcast  H  data  follow 
s  s 

the  trends  of  the  recorded  data  with  slight  phase  shifts  between  the  two 
data  sets.  If  the  phasing  were  adjusted,  the  differences  between  the 
data  would  be  about  ±0.5  ft. 

51.  Although  the  comparisons  indicate  good  agreement  between 

the  hindcast  and  measured  data,  extensive  comparisons  were  not  performed 
due  to  the  problem  of  the  gage-measured  wave  period.  Without  the  wave 
period  comparisons,  the  gage-measured  wave  data  are  not  appropriate  for 
use  in  verifying  the  SWWM. 

52.  In  order  to  provide  a  more  detailed  verification  of  the  SWWM, 
additional  comparisons  between  the  SWWM  and  measured  wave  data  are  shown 
in  Appendix  A.  The  data  are  obtained  from  a  similar  hindcast  study  in 
Saginaw  Bay,  Michigan.  The  methodology  of  the  SWWM  is  consistent  with 
the  Mississippi  Sound  study  and  Saginaw  Bay.  The  primary  purpose  for 
that  study  was  verification  of  the  SWWM  with  wave  gage  data.  Differ¬ 
ences  between  the  two  data  sets  are  ±0.5  ft  for  H  and  ±1.0  sec  for  T 

s 


*  N.  C.  Roberts.  1974.  "Wave  Height  Statistics,  Mississippi,"  Nash  C. 
Roberts,  Jr.,  Consultants  (typewritten  report). 
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(Station  4  location  is  found  in  Figure  8,  identified  as  WG  4) 


I 


(Garcia  and  Jensen,  in  preparation).  These  "error  bands"  associated 
with  the  hindcast  results  are  expected  to  apply  to  the  results  of  the 
Mississippi  Sound  study;  but  without  additional  field  data  (in  the 
Sound),  these  confidence  limits  cannot  be  taken  as  definitive. 


36 


PART  IV:  20-YEAR  HINDCAST  RESULTS 


53.  This  section  of  the  report  is  intended  only  to  serve  as  a 
general  description  of  simple  wave  characteristics  such  as  height,  pe¬ 
riod,  and  direction  of  wave  propagation.  More  detailed  analysis  of  one¬ 
dimensional  spectral  properties  and  interrelations  among  various  wave 
parameters  related  to  storm  characteristics  cannot  be  covered  in  the 
scope  of  this  report.* 

54.  The  20-year  hindcast  only  considers  a  single  population  of 
wave  conditions  defined  as  sea  (although  during  peiods  of  time,  the  non- 
dimensional  peak  frequency  is  less  than  0.13,  defined  as  swell  by  Hassel- 
mann  et  al.  (1976)).  Also,  all  wind  data  are  assumed  to  be  derived  from 
extratropical  storm  conditions;  therefore  wave  characteristics  calculated 
for  periods  of  tropical  storms  (Table  2)  must  be  regarded  as  approximate. 
Any  extremal  analysis  of  Hg  based  cn  these  data  will  probably  result 

in  an  underestimate  of  Hg  conditions.  All  water  depths  employed  in 
this  study  were  assumed  to  be  constant,  neglecting  tides  and  surges. 

55.  The  wave  parameters  Hg  (Equation  20),  T  (Equation  21), 
and  i()  given  at  each  station  every  3  hr  for  20  years  are  used  as  a 
basis  to  construct  all  products  found  in  Appendices  D-G.  The  value  of 
the  wave  direction  is  the  direction  from  which  the  waves  are  coming  (Fig¬ 
ure  14).  Azimuths  are  measured  clockwise  in  degrees  from  true  north  (0). 
As  previously  mentioned,  the  direction  of  wave  propagation  is  limited  to 
16  categories  (22.5-deg  increments  from  0  deg  to  337,5  deg). 

56.  Six  products  are  presented: 

a.  Seasonal  and  20-year  Percent  Occurrence  Tables. 

b.  Percent  Exceedance  Diagrams. 

c.  Height,  Period,  and  Direction  Histograms. 

d.  Mean  Hg  and  Largest  H  Tables. 

e.  "Over"  Duration  Tables. 

f.  "Under"  Duration  Tables. 


*  One-dimensional  frequency  spectra  and  Hg  ,  T  ,  and  4*  are  stored 
at  3-hr  intervals  for  20  years  for  each  station. 
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WAVE  DIRECTION 


180° 

Figure  14.  Diagram  indicating  wave  directions 


A  brief  description  of  each  product  is  given  and  sections  on  use  of  the 
products,  including  examples,  are  provided.  Each  product,  except  the 
duration  tables  (Appendices  F  and  G),  is  provided  separately  for  each 
station.  The  duration  tables  for  waves  over  and  under  specified  wave 
heights  are  presented  in  two  tabulations,  each  incorporating  all 
stations . 


Seasonal  and  20-Tear  Percent  Occurrence  Tables 


Description 


57.  Two  types  of 
directions  tables.  The 
currence  of  significant 


tables  are  printed:  azimuth  tables  and  all¬ 
azimuth  tables  (Figure  15)  give  the  percent  oc- 
waves  in  height  and  period  ranges  for  specif  ed 
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Figure  15.  Sample  percent  occurrence  tables  by  season  and  direction 


station,  season,"  ami  direction.  The  title  of  each  table  identifies  the 
station,  season,  angle  class,  and  water  depth.  The  wave  period  ranges 
are  in  1-sec  intervals  and  the  height  ranges  are  in  1/2-ft  increments. 
Values  in  the  azimuth  tables  represent  the  percent  of  the  season  that 
waves  occur  from  the  specified  angle  class  for  the  indicated  height  and 
period  ranges.  The  values  have  been  multiplied  by  1,000  to  allow  more 
accuracy  with  less  printing  space.  Summations  of  period  and  height 
ranges  are  provided  in  the  last  column  and  row  of  each  table.  The  sum¬ 
mations  are  also  multiplied  by  1,000.  The  last  line  in  each  directional 
table  contains  the  following  information  for  the  specified  angle  class 
and  season: 

a.  The  average  H 

s 

b.  The  largest  H 

s 

c.  Percent  of  waves  occurring  in  the  specified  season  from 
the  indicated  angle  class. 

58.  The  all-directions  table  for  each  season  (Figure  16)  is 
printed  after  the  337.5-deg  angle  class  table  for  the  specified  season. 
These  tables  give  the  percent  occurrence  of  significant  waves  within  the 
same  specified  height  and  period  ranges  coming  from  all  directions  for 
the  indicated  season  and  station.  Values  in  the  a  1 1 -d i rect i ons  tables 
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Figure  16.  Sample  percent  occurrence  table  for  one  season  and  all 

directions 


*  Season  1  is  December  through  February,  Season  2  is  March  through  May, 
Season  3  is  June  through  August,  and  Season  4  is  September  through 
November . 
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are  multiplied  by  100.  The  parameters  listed  in  the  last  line  of  the 

seasonal  tables  are  derived  from  all  directions  of  the  specified  season. 

The  number  of  cases  represents  the  number  of  wave  conditions  occurring 

in  the  indicated  season  for  the  total  20  years  analyzed.  All  calm  wave 

conditions  (H  <0.1  ft)  are  removed  from  these  tables  because  of  the 
s 

methodology  employed  in  the  SWWM.  This  constraint  will  not  alter  the 
data  since:  (a)  we  are  dealing  with  a  very  small  number  of  occurrences 
within  a  given  period  of  time,  and  (b)  the  magnitude  of  Hg  is  so  small 
that  whether  one  considers  a  wave  height  of  0.1  ft  or  0.0  ft  is  of  no 
real  consequence. 

59.  The  20-year  tables  (Figures  17  and  18)  also  included  in  this 
report  are  of  the  same  format  described  above.  These  tables  (angle 
class  tables  and  all-directions  tables)  are  based  on  the  full  20-year 
data  set  (or  the  cumulative  total  of  the  four  sets  of  seasonal  tables). 
Use  of  the  tables 

60.  The  tables  have  been  developed  to  produce  reliable  informa¬ 
tion  when  seasonal  or  20-year  intervals  are  employed. 

Examples 

61.  In  order  to  find  the  number  of  hours  that  waves  of  >1.0  ft  and 
<1.5  ft  and  3.0  to  3.9  sec  are  expected  to  occur  from  45  deg  at  Station  1 
in  Season  1  (December,  January,  February)  for  a  20-year  interval,  the 
value  read  in  the  table  for  the  specified  station,  season,  angle  class, 
height,  and  period  should  first  be  divided  by  1,000,  which  for  this  exam¬ 
ple  yields  0.006  percent  (Figure  15,  Appendices  D  and  E).  Then  0.006  is 
divided  by  100  to  give  the  probability,  and  multiplied  by  the  number  of 
hours  in  Season  1  for  the  20-year  interval  to  yield  the  number  of  hours 
that  the  specified  wave  is  expected  to  occur.  The  simple  conversion 
process  is: 


value  read  in  table  . 

1,000  T 


number  of 
100  x  hours  in 

time  interval 


number  of  hours 
=  specified  wave  is 
expected  to  occur 


(25) 


For  this  example: 


6 

1 ,000 


100  x  43,320  =  3  hr 
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STATION  1  20  YEARS  ANGLE  CLASS  (DEG  AZIMUTH  I  =  0. 

WATER  DEPTH  =  6.50  FEET 

PERCENT  OCCURRENCEIXIOOO >  OF  HEIGHT  AND  PERIOD  BT  DIRECTION 


HEIGHT(FEET) 


0.  -  0.49 

0.50  -  0.99 

}:S8 : 

2.00  -  2.49 

2.50  -  2.99 
3 . 00  - 

3 . 50  -  3.99 
4.00  -  4.49 

4.50  -  4.99 
5.00  -  GREATER 

TOTAL 


ITAL  0  2082 

AVERAGE  HS(FT)  =  0.47 


PERIOD! SECONDS ) 


3.0-  4.0- 

3.9  4. 9 


9.0- 

LONGER 


LARGEST  HS(FT)  =  1.64  ANGLE  CLASS  X 


STATION  1  20  YEARS  ANGLE  CLASS  (DEG  AZIMUTH)  =  22.5 

WATER  DEPTH  =  6.50  FEET 

PERCENT  OCCURRENCE! X1000 )  OF  HEIGHT  AND  PERIOD  BY  DIRECTION 


HEIGHT! FEET) 


_ 

0.49 

50 

- 

0.99 

oo 

- 

1.49 

50 

— 

1.99 

00 

- 

2.49 

*2 

— 

2.99 

00 

- 

3.49 

40 

- 

3.99 

on 

— 

4  .49 

40 

- 

4.99 

00 

- 

GREATER 

TOTAL 

AVERAGE  HS(FT) 


0  573 

=  0.45 


PERIOD! SECONDS) 

4.0-  5.0-  6.0- 

9  4.9  5.9  6.9 


LARGEST  HS(FT)  =  1.21 


9.0- 

LONGER 


0  0 
ANGLE  CLASS  X 


SimEPT&  =2Vf40PSFEE/NGU  CU5S  tDEG  «IWJTH 

PERCENT  OCCURRENCE ( X1000 )  OF  HEIGHT  AND  PERIOD  BY 


DIRECTION 


HEIGHT! FEET) 


0. 

- 

0.49 

0.50 

- 

0.99 

1.00 

— 

1.49 

1.50 

- 

1.99 

2.00 

- 

2.49 

2.50 

- 

2.99 

3.00 

- 

3.49 

3.50 

- 

3.99 

4.00 

— 

4.49 

4.50 

- 

4.99 

5.00 

GREATFP 

PERIOD! SECONDS) 


9.0- 

LONGER 


AVERAGE  HS(FT)  =  0.32  LARGEST  HS(FT)  =  1.05  ANGLE  CLASS  X  = 


STATION  1  20  YEARS  ANGLE  CLASS  (DEG  AZIMUTH)  =  67.5 

WATER  DEPTH  =  6.50  FEET 

PERCENT  OCCURRENCEIXIOOO )  OF  HEIGHT  AND  PERIOD  BY  DIRECTION 


HEIGHT! FEET) 


0.  -  0.49 

0.50  -  0.99 

IM : 

2.00  -  2.49 

2.50  -  2.99 
3.00  -  3.49 

3.50  -  3.99 
4.00  -  4.49 

4.50  -  4.99 
5.00  -  GREATER 

TOTAL 


AVERAGE  HS(FT) 


6  3880 

=  0.28 


PERIOD! SECONDS) 

4.0-  5.0-  6.0 

9  4.9  5.9  6 


2234  3  6 

LARGEST  HS(FT)  =  1.21 


9.0- 

LONGER 


ANGLE  CLASS  X  = 


Figure  17.  Sample  percent  occurrence  tables  for  all  seasons  (20  years) 

by  direction 


HATER  DEPtS^TsO  'FEE?0  YEARS  F*  AU  OIRECTIOf<S 
PERCENT  OCCURRENCE! X100 )  OF  HEIGHT  AND  PERIOD  FOR  ALL  DIRECTIONS 

HEIGHT! FEET)  PERIOD! SECONDS ) 


0. 

0.50 

1.00 

1.50 
2.00 

2.50 
3.00 

3.50 
4.00 

4.50 
5.00 


0.49 

0.99 

1.49 

1.99 

2.49 

2.99 

3.49 

3.99 

4.49 

4.99 
GREATER 


TOTAL 

AVE  HS(FT) 


0.0- 

0.9 

1.0- 

1.9 

2.0- 

2.9 

3.0- 

3.9 

4.0- 

4.9 

5.0-  6 

5.9 

0-  7.0- 

6.9  7.9 

8.0-  9 

8.9 

3267 

3031 

2368 

7 

27 

346 

182 

1 

77 

48 

39 

2 

82 

140  : 

58 

i  82 

10 

24 

6 

3267 

5406 

556 

166 

82 

198  92 

24 

=  0.45 

LARGEST  HS(FT)  = 

3.19 

TOTAL 

CASES  = 

58440 

TOTAL 


6325 

3013 

295 

40 

64 

39 

3 

0 

0 

0 
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Figure  18.  Sample  percent  occurrence  table  for  all  seasons 
(20  years)  and  all  directions 


The  following  tabulation  lists  the  approximate  number  of  hours  in  each 
season  for  1  year: 


Number  of  Hours  in 
_ Season _  Season  for  1  Year 

1  (D,  J,  F)  -2166 

2  (M,  A,  M)  -2208 

3  (J,  J,  A)  -2208 

4  (S,  0,  N)  -2184 


The  total  number  of  hours  in  the  time  interval  for  a  specified  season 
can  be  found  by  multiplying  the  number  of  hours  listed  in  the  above 
tabulation  by  the  number  of  years  required  for  a  given  problem. 

62.  The  all-directions  tables  can  be  used  in  a  similar  fashion. 
To  find  the  number  of  hours  waves  >1.0  ft  and  <1.5  ft  are  expected  to 
occur  within  a  year  for  Station  1  in  Season  1  for  all  directions  and 
periods,  divide  the  value  in  the  total  column  for  the  specified  Hg 
range  by  100,  which  yields  a  percent  of  0.36  (Figure  16).  Divide  0.36 
by  100  to  get  the  probability;  then  multiply  by  the  number  of  hours  in 
Season  1  for  1  year.  That  is: 


360  . 
100  T 


100  x  2,166  S  78  hr 
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63.  This  procedure  may  also  he  used  when  dealing  with  the  data  in 
the  20-year  angle  class  and  a l 1 -d i rect i 011s  tables.  The  number  of  hours 
in  the  time  interval  found  in  Equation  25  would  now  be  the  sum  for  all 
four  seasons  or  8,766  hours  per  year. 

Percent  Exceedance  Diagrams 

Description 

64.  The  percent  exceedance  diagrams  (Figure  19)  were  constructed 
for  all  months  in  the  20  years  of  hindcast  significant  wave-height  data. 
The  diagrams  have  a  log  scale  percent  axis  to  allow  greater  accuracy  for 
the  larger  wave-height  categories.  One  must  note  that  the  larger  wave- 
height  categories  will  be  affected  by  tropical  storm  conditions.  The 
calculations  performed  in  this  study  assumed  all  storm-generated  winds 
and  waves  were  derived  from  only  extratropical  conditions,  even  though 
tropical  storm  periods  were  included. 

Example 

65.  To  find  the  percent  of  time  that  a  3-ft  wave  is  exceeded  at 
Station  1,  locate  the  3-ft  intersect  with  the  curve  supplied  in  the  dia¬ 
gram  (Figure  19)  and  read  the  percent  exceeded  from  the  percent  axis 
(for  this  example,  ~0.025  percent). 

Height,  Period,  Direction  His tograms 


Description 

66.  Histograms  for  each  station  are  presented  for  seasonal  (Fig¬ 
ure  20)  and  full  20-year  data  (Figure  21).  The  histograms  show  the  per¬ 
cent  occurrence  of  significant  waves  as  a  function  of  the  three  basic 
wave  parameters  (Hg  ,  T  ,  and  tfj) . 

Use  of  the  diagrams 

67.  These  diagrams  are  intended  for  use  as  visual  aids  and  rela¬ 
tively  gross  estimations. 

Example 

68.  For  the  20-year  record  of  Station  1  wave  conditions, 
approximately  30  percent  of  the  waves  were  between  0.5  and  1.0  ft  in 
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•fa 


Figure  19.  Sample  percent  exceedance  diagram 
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Figure  20.  Sample  height,  period,  and  direction 
histograms  by  season 


Figure  21.  Sample  height,  period,  and  direction 
histograms  for  all  seasons  (20  years) 


height  (H^)  ,  about  6  percent  of  the  waves  had  peak  wave  periods  between 
3.0  and  4.0  sec  (T),  and  more  waves  came  from  0-deg  and  180-deg  angle 
classes  (4<)  than  any  other  category  (Figure  21). 

Mean  and  Largest  Hg  Table 

Description 

69.  Two  tables  which  summarize  the  mean  and  largest  for  each 

month  and  year  are  provided  for  each  station  (Figure  22).  The  mean  table 
also  provides  a  mean  monthly  value  and  a  mean  yearly  value  for  .  The 

largest  table  provides  the  largest  hindcast  for  each  year  and 

month  as  well  as  the  largest  hindcast  for  the  specified  station. 

Use  of  the  tab 1 e s 

70.  The  tables  can  be  used  as  a  quick  reference  in  determining 

gross  estimates  of  the  wave  climate  of  an  area.  Due  to  extreme  varia¬ 
tions  in  wave  heights,  the  mean  Hg  value  is  of  little  use  beyond  gross 
estimates.  The  largest  value  provides  an  idea  of  what  extreme  sig¬ 

nificant  wave  heights  might  occur;  however,  other  data  products  pre¬ 
sented  in  this  report  can  assist  in  determining  how  often  to  expect  ex¬ 
treme  values  and  how  long  they  might  occur. 

Example 

71.  To  determine  the  mean  Hs  at  Station  1  for  January  1956, 
simply  read  the  value  in  the  specified  column  and  row  (Figure  22).  The 
mean  Hs  for  1956  is  given  in  the  mean  column  opposite  1956.  The  mean 

for  all  January's  is  given  in  the  mean  row  under  January  for  this 
example : 


a . 

The 

mean 

H 

s 

for  January  1956  = 

0.5  ft. 

b. 

The 

mean 

H 

s 

for  1956  =  0.4  ft. 

c . 

The 

mean 

H 

s 

for  all  January's 

=  0.5  ft 

The  largest  Hs  table  can  be  read  in  a  similar  fashion,  and  by  scanning 
the  columns  and  rows,  additional  information  can  be  determined  (Fig¬ 
ure  22). 


a . 

The 

largest 

H 

s 

for  January  1956  =  2.5  ft. 

b. 

The 

largest 

H 

s 

for  1956  =  2.6  ft. 

c . 

The 

largest 

H 

s 

for  all  January's  =  3.2  ft 
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MEAN  HS(FEET)  BY  MONTH  AND  YEAR 
STATION  1 
MONTH 


YEAR 

JAN 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

OCT 

NOV 

DEC 

MEAN 

1956 

0.5 

0.4 

0.4 

0.4 

0.3 

0.4 

0.4 

0.3 

0.3 

0.3 

0.4 

0.4 

0.4 

1957 

0.4 

0.4 

0.4 

0.4 

0.4 

0 . 3 

0.4 

0.4 

0.3 

0.4 

0.5 

0.4 

0.4 

1958 

0.5 

0.5 

0.5 

0.6 

0.4 

0.4 

0.4 

0.4 

0.4 

0.5 

0.5 

0.6 

0.5 

1959 

0.5 

0.5 

0.6 

0.5 

0.4 

0.5 

0.4 

0.5 

0.5 

0.5 

0.6 

0.6 

0.5 

1960 

0.6 

0.7 

0.7 

0.5 

0.5 

0.4 

0.5 

0.3 

0.4 

0.4 

0.5 

0.5 

0.5 

1961 

0.6 

0.7 

0.6 

0.7 

0.5 

0.4 

0.4 

0.4 

0.4 

0.5 

0.6 

0.6 

0.5 

1962 

0.7 

0.7 

0.7 

0.5 

0.5 

0.4 

0.4 

0.3 

0.3 

0.3 

0.4 

0.4 

0.5 

1963 

0.4 

0.5 

0.4 

0.5 

0.5 

0.5 

0.4 

0.4 

0.4 

0.3 

0.5 

0.6 

0.5 

1964 

0.6 

0.7 

0.6 

0.6 

0.4 

0.4 

0.5 

0.4 

0.4 

0.6 

0.5 

0.6 

0.5 

1965 

0.6 

0.7 

0.7 

0.6 

0.5 

0.5 

0.4 

0.4 

0.4 

0.5 

0.4 

0.4 

0.5 

1966 

0.5 

0.6 

0.6 

0.7 

0.5 

0.4 

0.4 

0.3 

0.3 

0.3 

0.4 

0.5 

0.5 

1967 

0.4 

0.5 

0.5 

0.5 

0.5 

0.4 

0.5 

0.3 

0.3 

0.4 

0.5 

0.5 

0.4 

1968 

0.4 

0.5 

0.5 

0.4 

0.4 

0.4 

0.3 

0.4 

0.3 

0.3 

0.5 

0.5 

0.4 

1969 

0.5 

0.5 

0.6 

0.5 

0.4 

0.4 

0.4 

0.3 

0.3 

0.4 

0.4 

0.5 

0.4 

1970 

0.4 

0.5 

0.5 

0.5 

0.4 

0.5 

0.4 

0.4 

0.4 

0.5 

0.6 

0.6 

0.5 

1971 

0.4 

0.6 

0.7 

0.6 

0.5 

0.6 

0.5 

0.5 

0.3 

0.3 

0.4 

0.4 

0.5 

1972 

0.5 

0.4 

0.5 

0.5 

0.4 

0.6 

0.4 

0.4 

0.3 

0.4 

0.5 

0.4 

0.4 

1973 

0.5 

0  5 

0.5 

0.6 

0.6 

0.4 

0.3 

0 . 3 

0.4 

0.3 

0.4 

0.6 

0.5 

1974 

0.4 

0.5 

0.5 

0.5 

0.5 

0.4 

0.4 

0.3 

0.4 

0.3 

0.4 

0.5 

0.4 

1975 

0.4 

0.5 

0.5 

0.4 

0.3 

0.4 

0.4 

0.3 

0.3 

0.3 

0.4 

0.4 

0.4 

MEAN 

0.5 

0.6 

0.6 

0.5 

0.5 

0.4 

0.4 

0.4 

0.4 

0.4 

0.5 

0.5 

LARGEST  HS(FEET)  BY  MONTH  AND  YEAR 
STATION  1 
MONTH 


JAN 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

OCT 

NOV 

DEC 

YEAR 

1956 

2.5 

2.4 

2.5 

2.4 

1.1 

1.1 

2.5 

1.0 

1.2 

1.1 

1.4 

2.6 

1957 

2.4 

2.4 

2.7 

2.8 

2.4 

1.2 

2.4 

2.3 

1 . 0 

1.0 

2.4 

2.6 

1958 

2.8 

2.8 

2.6 

2.8 

2.4 

2.5 

2.7 

2.3 

1.0 

2.4 

1 . 7 

2.3 

1959 

2.5 

2.3 

2.6 

1.2 

1.5 

2.5 

1.1 

2.5 

2.4 

2.9 

2.4 

2.9 

1960 

3.2 

2.8 

3.1 

2.9 

2.4 

2.8 

2.4 

2.3 

1.8 

1.4 

1.2 

2.9 

1961 

3.1 

2.9 

2.9 

3.1 

2.6 

2.3 

2.4 

2.3 

1.3 

2.4 

2.9 

2.7 

1962 

3.0 

2.7 

3.0 

2.2 

2.3 

2.6 

2.3 

1.2 

1.2 

1.0 

2.4 

2.4 

1963 

1.7 

2.5 

2.6 

2.5 

2.6 

2.3 

2.4 

2.3 

1.5 

1.0 

3.0 

2.3 

1964 

3.0 

3.2 

3.0 

2.4 

2.3 

2.5 

2.7 

2.4 

1.0 

2.3 

2.4 

2.5 

1965 

3.0 

3.0 

2.9 

2.7 

2.6 

2.6 

2.7 

2.8 

1.2 

2.4 

2.6 

1.1 

1966 

2.6 

2.9 

2.8 

3.0 

1.2 

1.3 

2.3 

1.4 

1.4 

2.3 

1.5 

2.7 

1967 

2.4 

2.4 

2.5 

2.7 

2.4 

2.3 

2.4 

1.3 

1.2 

1.3 

2.3 

2.9 

1968 

2.9 

2.3 

2.8 

2.3 

2.5 

2.5 

1.1 

2.6 

1.0 

1 . 1 

2.5 

2.7 

1969 

2.3 

2.4 

2.9 

2.8 

2.3 

2.3 

2.3 

2.4 

1.2 

2.3 

2.4 

2.7 

1970 

2.3 

2.5 

2.7 

2.3 

1.2 

2.7 

2.8 

2.5 

2.3 

2.4 

2.7 

2.7 

1971 

2.7 

2.9 

3.0 

2.8 

2.7 

2.3 

2.3 

2.7 

1.3 

2.3 

2.4 

1.1 

1972 

2.4 

2.6 

2.9 

2.3 

1.2 

2.7 

2.4 

2.4 

2.5 

2.3 

2.6 

2.6 

1973 

2.5 

2.6 

2.7 

2.6 

2.5 

2.4 

1.2 

2.3 

1.4 

1.6 

2.8 

2.8 

1974 

2.4 

2.6 

2.5 

2.3 

2.4 

2.3 

2.5 

1.2 

1.2 

0.8 

2.6 

2.6 

1975 

2.3 

2.5 

2.3 

1.2 

1.1 

2.6 

1.4 

2.4 

1.5 

1.4 

1.4 

2.7 

LARGEST  HSCFEET)  FOR  STATION  1  =  3.2 

Figure  22.  Sample  mean  Hg  and  largest  tables 
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Duration  Tables 


Description 

72.  These  tables  contain  H  duration  information  for  all  23  Mis- 

s 

sissippi  Sound  stations.  The  values  given  in  the  tables  are  the  mean 
duration  in  hours  (x)  and  the  maximum  duration  in  hours  (mx).  For  the 
"over"  tables  (Figure  23  and  Appendix  F),  duration  is  defined  as  the 
length  of  time  waves  greater  than  a  given  height  persist  once  the  wave 
height  has  been  exceeded,  and  for  the  "under"  tables  (Figure  24  and  Ap¬ 
pendix  G) ,  duration  is  the  length  of  time  waves  less  than  a  given  height 
persist  once  the  waves  become  less  than  the  indicated  Hg  category. 

The  computations  used  data  from  all  months  of  the  20-year  hindcast  wave 
data . 

Use  of  the  tables 

73.  For  sequences  of  wave  heights,  discretized  sets  of  durations 
are  formed  of  wave  heights  above  or  below  a  specified  level  after  waves 
exceed  or  become  less  than  that  level.  Each  element  of  an  individual 
set  is  a  single  duration  event,  and  the  expected  duration  above  or  below 
a  specified  level  is  equal  to  the  mean  of  all  elements  in  the  set. 

Thus : 


n  . 
J 


i=l 


where 


subscript  j  =  an  incremented  wave-height  counter  for  wave  height 

h. 

J 

i  =  a  counter  for  duration  events 

n  =  number  of  times  that  waves  exceed  or  become  less 
than  a  particular  wave  height 

x.  .  =  duration  of  a  single  event  with  wave  heights  above 

^  or  below  level  H. 

J 

In  this  formulation  a  number  of  other  useful  duration  parameters  from 
the  20-year  sample  can  be  estimated. 
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Duration  of  Waves  over  a  Specified  Height 


Wave 

Height 

Stations 

1 

2 

3 

4 

5 

6 

7 

8 

Class 

ft 

x 

mx 

X 

mx 

X 

mx 

X 

mx 

X 

mx 

X 

mx 

X 

mx 

X 

mx 

>0.5 

9 

117 

14 

>1 

14 

19S 

14 

261 

16 

273 

14 

24  9 

12 

174 

9 

150 

>1.0 

5 

30 

tj 

5** 

7 

81 

7 

90 

7 

114 

7 

114 

6 

93 

6 

78 

>1.5 

4 

19 

5 

30 

S 

19 

6 

45 

5 

42 

5 

57 

5 

42 

5 

27 

>2.0 

4 

15 

15 

4 

18 

4 

21 

4 

18 

4 

10 

4 

18 

4 

15 

>2.5 

4 

15 

4 

15 

4 

6 

4 

19 

4 

15 

4 

15 

4 

15 

3 

9 

>3.0 

3 

6 

4 

15 

3 

3 

4 

15 

4 

15 

4 

15 

4 

15 

3 

3 

>3.5 

— 

— 

4 

15 

-- 

-- 

4 

15 

4 

15 

4 

15 

4 

15 

— 

-- 

.*4.0 

— 

— 

4 

12 

-- 

-- 

4 

15 

4 

15 

4 

15 

3 

6 

— 

— 

>4.5 

— 

— 

3 

3 

— 

— 

4 

12 

4 

9 

4 

15 

— 

— 

— 

-- 

>5.0 

— 

— 

— 

— 

-- 

— 

3 

6 

3 

3 

3 

3 

— 

— 

— 

— 

>5.5 

— 

— 

— 

— 

— 

— 

3 

3 

-- 

— 

— 

— 

-- 

— 

-- 

— 

>6.0 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

>6.5 

— 

— 

-- 

-- 

-- 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

>7.0 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

(cont inued) 


Note:  x  equals  the  mean  duration  in  hours  for  the  specified  wave-height  class;  mx  equals  the  maximum 

duration  in  hours  for  the  specified  wave-height  class. 

Figure  23.  Sample  "over"  duration  table 


Duration  of  Waves  Under  a  Specified  Height 


Wave 

Height 

Class 

ft 

Stations 

1 

■> 

3 

4 

5 

X 

mx 

X 

mx 

x 

mx 

> 

< 

mx 

X 

mx 

Calm 

3 

3 

3 

3 

3 

3 

3 

6 

3 

3 

<0.5 

19 

273 

12 

135 

11 

165 

a 

180 

9 

96 

<1.0 

105 

1,365 

53 

804 

37 

804 

37 

633 

24 

492 

1.5 

246 

3,447 

185 

2,727 

165 

3,111 

121 

2,277 

91 

981 

<2.0 

299 

3,447 

273 

3,099 

932 

8,682 

236 

2,727 

186 

2,244 

<2.5 

1,113 

6,969 

293 

3,111 

6,64  9 

40,185 

280 

3,111 

212 

3,084 

<3.0 

8,178 

35,484 

294 

3,111 

53,218 

114,498 

293 

3,111 

214 

3,084 

<3.5 

— 

— 

84  3 

6,213 

-- 

— 

414 

6,165 

577 

t> .  1  f'  5 

<4.0 

— 

— 

2,809 

14,721 

— 

— 

1 

,157 

8,241 

1,619 

17,769 

<4.5 

— 

— 

11,828 

35,484 

— 

— 

2 

,809 

14,721 

3,639 

33,180 

<5.0 

— 

— 

— 

— 

— 

— 

7 

,384 

42,978 

11,828 

3>  ,484 

<5.5 

<6.0 

— 

— 

__ 

_ 

— 

35 

,490 

35.742 

__ 

<6.5 

<7.0 

:: 

_ 

__ 

_ 

:: 

_ 

__ 

__ 

(continued) 


Note:  x  equals  the  mean  duration  in  hours  for  the  specified  wave-height  class;  mx  equals  the  maximum 
duration  in  hours  for  the  specified  wave-height  class. 

Figure  24.  Sample  "under"  duration  table 


(1 

J 
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a.  x.  -  the  maximum  duration  of  waves  above  or  below  H.  . 

jmax  j 

b.  x.  .  -  the  minimum  duration  of  waves  greater  than  or 

jmin 

less  than  H. 

J 

c.  p(x^)  ”  t*le  probability  of  occurrence  of  a  particular 
duration  of  wave  heights  greater  than  or  less  than  H.  . 

d.  F(x^)  "  the  cumulative  probability  of  occurrence  of  dura¬ 

tions  less  than  or  greater  than  x  for  a  particular  H. 
category.  ^ 

e.  F'(Xj)  -  the  cumulative  probability  that  a  duration  event 
in  category  1L  will  exceed  or  be  less  than  length  x  . 

74.  In  particular,  the  parameter  F'(x.)  plays  a  significant 
role  in  the  assessment  of  design  criteria  when  such  criteria  are  based 
on  expected  damage  as  a  function  of  duration  of  given  wave  conditions. 

For  example,  if  the  wave  height  associated  with  a  20-year  return  period 
has  an  expected  duration  of  3  hr,  and  it  takes  a  6-hr  duration  to  sig¬ 
nificantly  damage  a  structure,  it  becomes  relevant  to  determine  the  prob¬ 
ability  of  a  6-hr  (or  greater)  duration  at  this  wave  height.  In  this 
sense,  the  design  wave  height  would  become  a  function  of  both  wave  height 
and  duration,  and  optimization  of  design  would  require  the  consideration 
of  both  F(x^)  and  F 

Example 

75.  To  find  the  mean  and  maximum  duration  for  waves  greater  than 
2.0  ft  and  less  than  2.5  ft  for  Station  1,  look  up  Station  1  in  the 
"over"  duration  table  and  read  the  values  in  the  >2.0  ft  row.  For  this 
example,  4  hr  is  the  mean  and  15  hr  is  the  maximum  (Figure  23). 


*  Since  this  is  calculated  only  when  is  exceeded,  this  number  is 

not  zero;  however,  for  most  intermediate  to  high  wave  heights,  the 

value  of  x.  .  is  equal  to  the  duration  associated  with  one  sample, 
jmin 

Only  for  small  waves  does  this  parameter  have  significance. 
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PART  V:  DISCUSSION  OF  RESULTS 


76.  This  report  presents  the  three  major  phases  of  the  Mississipp 
Sound  Wave-Hindcast  Study: 

a.  Estimation  of  long-term  wind  conditions  overwater. 

b.  Development  of  a  parametric  shallow-water  wave  model. 

c.  Calculation  of  long-term  wave  statistics  for  the  study 
area . 

Each  phase  of  the  study  was  treated  in  detail,  since  errors  in  any  phase 
could  alter  the  results  in  subsequent  phases  dramatically.  The  estima¬ 
tion  of  wind  conditions  in  Mississippi  Sound  from  overland  wind  condi¬ 
tions  employed  techniques  that  have  been  verified  through  extensive 
field  studies  (Resio  and  Vincent  1976).  The  only  assumption  affecting 
the  wind  results  was:  for  each  time  interval,  wind  speeds  and  direc¬ 
tions  were  uniform  over  the  Sound.  This  assumption  was  verified  in 
PART  II  of  this  report. 

77.  The  basis  for  the  development  of  the  Shallow-Water  Wave  Model 
(SWWM)  was  to  adequately  describe  the  physical  process  involved  in  wave 
growth  and  finite  depth  wave  transformations  while  maximizing  computa¬ 
tional  efficiency.  Detailed  comparisons  of  the  SWWM  results  with  mea¬ 
sured  and  visual  observations  in  Mississippi  Sound  are  necessary  for 
verification.  Unfortunately,  only  limited  wave  data  were  available  in 
Mississippi  Sound  to  perform  comparisons;  and  verification  of  the  mod¬ 
eling  technique  is  not  possible  with  such  limited  information.  Results 
from  another  study  employing  the  SWWM  (Appendix  A  and  Garcia  and  Jensen, 
in  preparation)  were  presented  and  clearly  showed  that  the  hindcast  data 
compare  favorably  with  the  measured  wave  data.  The  physical  processes 
and  geometry  associated  with  both  studies  are  very  similar  (excluding 
tropical  storms).  Therefore,  it  is  reasonable  to  assume  the  SWWM  ade¬ 
quately  models  finite  water  depth  wave  growth  and  transformation  pro¬ 
cesses  and  that  application  to  Mississippi  Sound  is  valid. 
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APPENDIX  A:  VERIFICATION  OF  THE  SHALLOW-WATER  WAVE  MODEL 


1.  The  shallow-water  wave  model  (SWWM)  was  first  applied  to  a 
project  in  Saginaw  Bay,  Michigan,  because  wave  data  were  already  avail¬ 
able  for  calibration  and  verification  in  that  location.  The  bathymetry 
(except  for  navigational  channels)  and  basin  configuration  of  Mississipp 
Sound  are  somewhat  similar  to  that  in  Saginaw  Bay  (Figure  Al).  Both 
basins  are  virtually  cut  off  from  adjacent  expansive  deepwater  areas  by 
shoals  and  islands.  The  assumptions  and  processes  involved  in  the  SWWM 
were  the  same  for  both  the  Saginaw  Bay  and  Mississippi  Sound  studies. 

Four  storm  conditions  were  selected  from  the  Saginaw  Bay  study  to  show 
the  relative  accuracy  of  the  modeling  technique.  The  scarcity  of  wave 
measurements  in  the  Mississippi  Sound  region  have  resulted  in  the  use  of 
data  from  Saginaw  Bay  (having  physical  and  geographical  structures  sim¬ 
ilar  to  the  Sound)  to  verify  the  SWWM.  Although  tb  s  is  an  indirect 
verification  of  the  model  for  its  use  on  Mississippi  Sound.,  it  is  the 
only  procedure  available  to  provide  a  reasonable  assurance  that  the 
model  is  valid  for  this  application. 

2.  Measured  wave  and  wind  data  from  the  9-11  May  1981  storm  were 
used  to  evaluate  and  verify  the  modeling  technique  in  Saginaw  Bay.  The 
Corps  anemometer  was  located  approximately  12  miles  southwest  of  Sta¬ 
tion  1.  This  storm  (Figure  A2)  produced  the  largest  waves  recorded 
since  initiation  of  that  measurement  program.  These  data  were  used  to 
calibrate  the  SWWM  within  Saginaw  Bay.  A  set  of  fetch  lengths  and  water 
depths  dependent  on  discrete  wind  direction  angles  were  computed  for 
each  station  (see  Figure  Al  for  the  station  locations).  The  methodology 
used  in  the  wave  model  was  the  same  for  both  stations.  Comparisons  be¬ 
tween  the  computed  and  measured  Hg  (significant  wave  height  defined  in 
Equation  21)  at  Stations  1  and  2  are  shown  in  Figures  A3  and  A4.  Compar¬ 
isons  between  the  computed  and  measured  T  (peak  period,  defined  in 
Equation  22),  are  shown  in  Figures  A5  and  A6.  The  SWWM  results  compare 
very  favorably  with  the  measured  Hg  data  for  both  stations;  and  the  T 
results,  over  an  average,  also  compare  favorably  with  the  measured  data. 
One  limitation  to  the  SWWM  is  that  it  assumes  discrete  frequency  bands 
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of  greater  width  than  the  results  obtained  by  the  analysis  of  the  gage 
records.  Large  variations  in  the  peak  period  are  a  function  of  time  and 
cannot  be  resolved  given  near-constant  wind  conditions. 

3.  Three  additional  comparisons  are  performed  and  shown  in  Fig¬ 
ures  A7-A13.  The  first  comparison  BCT  1  (Blind  Comparison  Test)  re¬ 
quired  a  secondary  source  of  wind  data  because  the  anemometer  deployed 
near  Saginaw  Bay  was  not  operable  during  the  outline  test  period  (5  and 
6  May  1981).  The  secondary  anemometer  was  at  Wurtsmith  Air  Force  Base 
located  near  Oscoda,  Michigan  (approximately  55  miles  northeast  of  Sta¬ 
tion  1).  The  average  wind  condition  for  BCT  1  is  approximately  15  mph, 
with  maximum  winds  of  25  mph  (Figure  A7).  The  wind  direction  is  shown 
to  be  highly  variable  during  this  particular  storm  (Figure  A7).  The 
results  displayed  in  Figures  A8  (Hg)  and  A9  (T)  show  how  capable  the 
SWWM  is  in  computing  accurate  wave  conditions  under  variable  wind  condi¬ 
tions.  The  next  set  of  comparisons  (BCT  2),  for  low  wind  speed  condi¬ 
tions  and  variable  wind  directions  (Figure  A10),  shows  that  the  SWWM  will 
predict  the  largest  recorded  Hs  (~2.5  ft)  in  the  event  and  approximate 
the  associated  T  (~4.0  sec,  Figure  All).  There  is  an  apparent  diver¬ 
gence  in  the  two  data  sets  after  the  maximum  condition  is  attained,  but 

the  differences  are,  at  its  maximum,  1.0  ft  for  H  and  0.6  sec  for 

s 

T  .  The  variability  in  the  wind  direction  under  low  wind  speeds  could 

cause  these  differences.  The  last  set  of  comparisons  BCT  3  (Figures  A12 

and  A13)  is  for  near-constant  wind  directions  and  moderate  wind  speeds. 

Figure  A12  displays  two  sets  of  wind  information,  the  Corps  anemometer 

(located  near  Bay  City,  approximately  10  miles  southwest  of  Station  1) 

and  Wurtsmith  Air  Force  Base.  One  can  see  a  major  phase  difference  in 

the  two  data  sets  of  approximately  14  hr.  Using  the  Corps  anemometer 

results  for  wind  condition  data  in  the  model,  we  find  that  the  hindcast 

wave  results  (Figure  A13)  for  both  Hg  and  T  lag  the  measured  results 

by  the  same  14  hr.  The  storm  event  is  correctly  simulated  by  the  SWWM, 

where  maximum  differences  fall  within  ±0.5  ft  for  H  and  ±0.5  sec  for 

s 

T  ,  if  the  hindcast  data  are  shifted  14  hr. 

4.  The  data  presented  in  this  discussion  are  only  examples  of 
many  comparisons  made  between  measured  and  hindcast  results.  The  only 
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changes  made  in  these  comparisons  were  the  input  wind  conditions, 
adjustments  in  coefficients  or  mechanisms  were  made  in  the  SWWM. 
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Figure  A4.  Comparison  between  measured  and  hindcast  significant 
wave-height  data  for  Station  2 
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Figure  A5.  Comparison  between  measured  and  hindcast  peak  wave 
period  data  for  Station  1 
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Figure  A6.  Comparison  between  measured  and  hindcast  peak  wave 
period  data  for  Station  2 
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Figure  A7.  Wind  information  data  for  BCT  1 
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igure  A8.  Comparison  between  measured  and  hindcast  significant 
wave-height  data  for  Station  1  for  BCT  1 
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Figure  A9.  Comparison  between  measured  and  hindcast  peak  wave 
period  data  for  Station  1  for  BCT  1 
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Figure  A10.  Wind  information  data  for  BCT  2 
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Figure  A12.  Wind  information  data  for  BCT  3 
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Figure  A13.  Comparisons  between  measured  and  hindcast  significant 
wave  height  and  peak  period  results  for  Station  1 ,  BCT  3 
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APPENDIX  B:  NOTATION 


E 

E(f) 

E(fj) 

E 


ff 

f 

c 
f . 
Jt 
f 

m 
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t 

F(x . ) 


F*(x  ) 


g 

h 

h. 

1 


m 

ti 

s 

i 


subscript  j 


k. 

1 


L. 

1 

mx 


Final  total  energy  at  F^ 

Continuum  energy  density 
Discrete  energy  density 

Total  energy  resulting  from  a  wind  speed 
Residual  energy 

Original  total  energy  at  F^_} 

Resultant  discrete  energy  density  spectrum  at  a  given 
station 

Fraction  factor 

Lower  frequency  bounding  the  total  energy 
Discrete  frequency  band 
Nondiraensional  peak  frequency 
Fetch  length 
Total  fetch  length 

The  cumulative  probability  of  occurrence  of  durations 
less  than  or  greater  than  x  for  a  particular  H. 
category  J 

The  cumulative  probability  that  a  duration  event  in  cate¬ 
gory  Hj  will  exceed  or  be  less  than  length  x 

Acceleration  due  to  gravity 

Mean  water  depth  along  F 

Water  depth  at 

Water  depth  at  a  given  station 
Wave  height 

Maximum  wave  condition 
Significant  wave  height 
Increment  counter 
Incremented  wave-height  counter 
Wave  number 

Nonvarying  parameter  (and  constants) 

Shoaling  coefficient 
Wavelength 
Maximum  duration 
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n . 
J 

p(xj) 

Q 

t 

t  . 
min 

T 

T 
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U 
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X  .  . 

1J 

X  . 

jmax 

x .  . 

jmin 

/■v 

X 

Of 

y 

AF. 

1 
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4>(wh) 

♦f 
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Number  of  times  waves  exceed  or  become  less  than  a  par¬ 
ticular  wave  height 

The  probability  of  occurrence  of  a  particular  duration  of 
wave  heights  greater  than  or  less  than  fT 

Dependent  parameter 

Total  elapsed  time  since  wind  began  to  blow 

Minimum  duration  condition 

Peak  wave  period 

Significant  wave  period 

Wind  speed 

Mean  duration 

Duration  of  a  single  event  with  wave  heights  above  or 

below  level  H. 

J 

The  maximum  duration  of  waves  above  or  below  H. 

J 

The  minimum  duration  of  waves  greater  than  or  less  than 

H. 

J 

Nondimens ional  fetch  length 
Phillips'  equilibrium  constant 
Weighting  function 

Distance  of  wave  travel  within  discrete  fetch  length 

Independent  parameter 

Wind  direction 

Nondimens ional  function 

Friction  function 

Direction  of  wave  propagation 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  I.ihrary  of  Congress  MARC  format  is  reproduced 
below. 
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